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We examine the entanglement induced by an angular momentum coupling between two harmonic 
systems. The Hamiltonian corresponds to that of a charged particle in a uniform magnetic field 
in an anisotropic quadratic potential, or equivalently, to that of a particle in a rotating quadratic 
potential. We analyze both the vacuum and thermal entanglement, obtaining analytic expressions 
for the entanglement entropy and negativity through the gaussian state formalism. It is shown that 
vacuum entanglement diverges at the edges of the dynamically stable sectors, increasing with the 
angular momentum and saturating for strong fields, whereas at finite temperature, entanglement is 
non-zero just within a finite field or frequency window and no longer diverges. Moreover, the limit 
temperature for entanglement is finite in the whole stable domain. The thermal behavior of the 
gaussian quantum discord and its difference with the negativity is also discussed. 



PACS numbers: 03.65.Ud,03.67.Mn,05.30. Jp 



I. INTRODUCTION 

The investigation of quantum entanglement and quan- 
tum correlations in distinct physical systems is of great 
interest for both quantum information and many-body 
physics [ll-ll]- While the evaluation of entanglement 
in systems with a high dimensional Hilbert space is in 
general a difficult problem, boson systems described by 
quadratic Hamiltonians in the basic boson operators of- 
fer the invaluable advantage of admitting an exact eval- 
uation of entanglement measures in both the ground and 
thermal state, through the gaussian state formalism [1,0- 
[T^ . The latter allows to express the entanglement en- 
tropy [m and negativity [ij, [ll| of any bipartition of 
a gaussian state in terms of the symplectic eigenvalues 
of covariance matrices of the basic operators. Moreover, 
the positive partial transpose (PPT) separability crite- 
rion [l3,[i3| is both necessary and sufficient for two-mode 
mixed gaussian states [1] (and also (l,n — 1) bipartitions 
of n modes gaussian states 0), turning the negativity 
into a rigorous entanglement indicator for these systems. 
Let us also remark that there is presently a great interest 
in continuous variable based quantum information [Tsj . 
where gaussian states constitute the basic element. 

In addition, an approximate yet analytic evaluation of 
the quantum discord [l9 . [20l in two- mode gaussian states 



was recently achieved [21|, [22|, by restricting the local 
measurement that determines this quantity to a gaus- 
sian measurement [l^ . Quantum discord is a measure of 
quantum correlations which coincides with the entangle- 
ment entropy in pure states but differs essentially from 
entanglement in mixed states, where it can be non-zero 
even if the state is separable, i.e., with no entanglement. 
The current interest in the quantum discord was trig- 
gered by its presence [23| in certain mixed state based 
quantum computation schemes which provide exponen- 
tial speedup over classical ones [l^], yet exhibiting no 
entanglement [25] . Important properties of states with 
non-zero discord were recently unveiled 26- 28)] . 

The aim of this work is to examine, using the gaussian 
state formalism, the entanglement and quantum corre- 



lations between two harmonic modes generated by an 
angular momentum coupling. Such system arises, for in- 
stance, when considering a charged particle in a uniform 
magnetic field in an anisotropic quadratic potential, or 
also a particle in a rotating anisotropic harmonic trap 
[29I - I32I . The model has then been employed in several 
areaSj including the description of deformed rotating nu- 
clei [31i, 32] , anisotropic quantum dots in a magnetic field 
[33} , and fast rotating Bose-Einstein condensates [s^ - fSTj 
in the lowest Landau level approximation [ssHioj . Con- 
taining just quadratic couplings in the associated boson 
operators, the different terms in the Hamiltonian may 
in principle be also simulated by standard optical means 
[3,|4l|. For a general quadratic potential, the model ex- 
hibits a complex dynamical phase diagram [4^ , present- 
ing distinct types of stable and unstable domains and 
admitting the possibility of stabilizing an initially un- 
stable system by increasing the field or frequency i43]. 
The model provides then an interesting and physically 
relevant scenario for analyzing the behavior of mode en- 
tanglement in different regimes and near the onset of 
different types of instabilities, with the advantage of al- 
lowing an exact analytic evaluation of entanglement and 
quantum correlation measures at both zero and finite 
temperature. In addition, the present results indicate 
that mode entanglement can be easily controlled in this 
systems by modifying the field or frequency, suggesting 
a potential for quantum information applications. Let 
us finally mention that the dynamics of entanglement in 
other two-mode systems were examined in ^43i-45.] . 

In sec. In] we describe the model and derive the ana- 
lytic expressions for the vacuum entanglement entropy 
and the thermal negativity. The basic features of the 
quantum discord are also discussed. The detailed behav- 
ior of entanglement with the relevant control parameters 
is then analyzed in sec. IIIIl where we show that while 
vacuum entanglement diverges at the edges of stable sec- 
tors, being correlated with the angular momentum, at 
finite temperature entanglement is finite, and non-zero 
just within a finite field window and below a finite limit 
temperature. A comparison between the thermal behav- 
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ior of the negativity and that of the gaussian quantum 
discord is finally made, which indicates a quite different 
thermal response of these two quantities, with the discord 
vanishing only asymptotically for T — > oo. Conclusions 
are finally drawn in IIVI 



B. Diagonalization and stability 



If the parameter 



A 



fc;)V4 



(6) 



II. FORMALISM 



Model Hamiltonian 



We consider a system described by the Hamiltonian 



H 



1 



KQi) + i^{p; + KQy) - ^(Q-Py - QyP^^) 



which represents two harmonic modes coupled by an 
angular momentum term. Here Q^i, stand for di- 
mensionless coordinates and momenta {[Qfj,,P,j] — iS^^, 
[<5/i,Qi/] — [Pti,Pu] = 0). Eq. (P) arises, for instance, in 
the description of a particle of charge e and mass m in 
a general quadratic potential subject to a uniform mag- 
netic field, parallel to a principal axis of the potential. 
Denoting this axis as z, such Hamiltonian reads 



n = 



{P~eA/cf 1 



2m 



KyQl 



K.Ql) (2) 



KQl + K'Ql-nc,]+n,, (3) 



e|H| 



the 



where A = x Q is the vector potential, ft 

cyclotron frequency, £3 = QxVy - QyVx, '^(■z = 2(1^ + 
KzQl) and K'^ = Kf, + mfl^/A. Eq. © is also identical 
with the intrinsic Hamiltonian of a particle of mass m 
in a quadratic potential of constants K'^ rotating around 
the z axis with frequency /2 [sil, HI] . 

The replacement Vf^ — P^-s/HtoHo, = 
Qfj,/ \/ milo/h in ([3|), with a reference frequency, leads 
toH — MIqH + "Hz, with H given by (P) and 



k' 



n 



(4) 



We note that in terms of the boson operators 6^ = 

{Qii+iPfj.)/ V2, the scaled angular momentum Lz = Cz/h 
in (P is 



Lz = QxPy - QyPx = -iip^by - blhx) , 



(5) 



and can then be simulated by standard linear optics, al- 
though for such bosons the first two terms in ITJ become 



Y.,^x.y9t{blb, + \)+g-{bl + bl ),with5, 



± _ 



k'±i 



, and 



J fi=x ,y d fj, \" fi" ' 21 ' illi\"fJ, ' "11 " '^"^'^ U fj. 2 

require non-linear means. If K'^ > 0, we can set k'^ = 1, 
i.e., g~ = 0, by adequately fixing fio in (|3]), but g~ re- 
mains non-zero in the relevant anisotropic case ky 7^ fc^, 
where [H,Lz] 7^ 0. The change to normal x,y bosons, 
such that H = J2fj,=x.y \/K^(^i^f^ + ^) ^ ^^z' ^^^^ ^^^'^ 



instead to an additional term oc (1 — y -jp-)ibxby "x"y/ 
in L,. 



blb\) 



is non-zero, the canonical transformation 

where 7 = —^7—^, V — ^^'^ labels {x,y) are 

now identified with — ), allows to write Eq. ([T|) as [i^l 

^= E (8) 



where a± = 1 - ^(7 T w) and /?± = ^(7 ± w). If 
A = and w 7^ 0, a separable representation of the 
form ([5]) in terms of canonical variables {\Q\,, P^] = Sfiv, 
[P'^^Pl] — = 0) is not feasible. Such a possi- 

bihty can arise in the repulsive case fc^ < 0, when the 
4x4 matrix representing the quadratic form (P is not 
diagonalizable with the symplectic metric and leads to 
non-trivial Jordan forms [42|. 

For general real values of k'^ in (P, the coefficients 
a^, 13^ in dp can be positive, zero, negative, and even 
complex [431 . We will here consider those cases where 
Eq. ([S]) can be further written as 



|A±| = \Ja±l3± = 



K+K 



with real and 6'^ 




(9) 

-u;2±A, (10) 



P' standard 



such 



fj,7 ^ Vi 0), 

In these cases the 



bosons ([6'^, b' 

that H exhibits a discrete spectrum 
matrix representing (P is diagonalizable with the sym- 
plectic metric, with real symplectic eigenvalues j42j . 

At fixed kfj, in (charged particle in a magnetic field) , 
Eq. ([9]) is vahd in the following domains [42*1 (Fig. P: 

(A) kx > 0, ky > 0, where a± > 0, /3± > and A± > 
0. This is the standard case of an attractive quadratic 
potential, where H is positive definite and hence fully 
stable. 

(B) kx < 0, ky < 0, and 



\uj\ > UJc 



(11) 



where a+ > 0, /3+ > but «_ < 0, /?_ < 0, implying 
A+ > but A- < in ([5]). This is the case of a repulsive 
quadratic potential, where H becomes equivalent to a 
standard plus an inverted oscillator if |w| > c^c'- It has 
no minimum energy, but is dynamically stable, as the 
motion remains bounded (4^ . The dynamics around a 
quadratic potential maximum can then be stabilized by 
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FIG. 1. Left: Sectors with discrete spectrum at fixed fc^ 
(particle in a magnetic field plus quadratic potential) in the 
anisotropy-fi-equency plane ky /\kx\,Lo /uo, where ljq = y^ffcTf- 
A is the positive definite sector, corresponding to an attractive 
potential A;^ > for = x,y, whereas B the non-positive 
sector, corresponding to a repulsive potential fc^ < for /i = 
x,y with |tj| > Wc (Eq. pifl V Right: Same sectors at fixed 
k'fj, > Q (particle in a rotating quadratic potential) in the 
k'y / k'n^ , u / uiQ plane, where u'q = and we have set k'^ > 0. 
A is the positive definite sector (Eq. (I12p l. B\ the non-positive 
sector with fc^^ > (Eq. ([13])) and B2 that with fc^ < (Eqs. 
(|14|| - (|15|) ). Quantities plotted in all figures are dimensionless. 



a sufficiently strong field. 

(C) kx ^ ky = and lu ^ {Landau case), where A = 
2a;^, a+ = 1, /3+ = iuj'^ whereas «_ = /3- =0, leading to 
A+ = 2\uj\ and A_ = 0. Eq. (jH]) becomes a standard plus a 
vanishing oscillator. Here the final choice of b'_ , b'l_ is not 
fixed, as A_ = 0. We will set = ^^J(3'_ + iPL/V^, 
according to the ^ limit of the isotropic case kx = 
ky = k, where A± = ± cj, A = 2ijj\fk! , 7 = \fk! = 
1/7?, a^, = iA^/VF, = 2\^,^ and k' ^ k + . 
k' 



At fixed fc^ in ([T]) (rotating potential), the previous 
sectors are seen quite differently. Sector (A) corresponds 
to fc^ > 0, fc^ > and 



(12) 



|a;| < = Mm[^yk 



indicating a maximum allowable frequency in an attrac- 
tive rotating quadratic potential (right panel in Fig. [1} . 
Sector (B) corresponds to 



\uj\ > = Max[ 



(13) 



if k'^ > 0, k'y > 0. Thus, as the frequency \uj\ is increased 
above uj[.i a finite instability interval cu'^i < |w| < 
arises in the anisotropic case k'^ ^ k'y, although dynami- 
cal stability is again recovered for > addition, 
sector (B) also corresponds here to A:^ > and A:^ < 
(or viceversa), provided [43 | 

(14) 

-3fc; <k'y< -45) 



-k'^<k'y<0, 



Uj\> 
< |w| < UJ',^ 



Hence, a quadratic potential repulsive in one of the axes 
can be stabilized by increasing the frequency above 
if — 3fc^ < k'y < although stability holds just within a 
finite interval if — 3fc^ < k' < ~-k'^. Finally, the Landau 



C. Covariance matrix 

Both the vacuum |0') of the primed bosons in ([9]), 
and the thermal state 



p = Z-^ exp[-/3iJ] , 



(16) 



well defined in the stable region (A) (with /3 — 1/T > 
and Z = Tt exp[-/37J] = J^^^^ 2sinh(^^A^/2) )' are gaus- 
sian states [1, 04i3- Any expectation value, and in par- 
ticular the entanglement between the x and y modes 
in these states, will then be completely determined by 
the elements of the basic covariance matrix of the op- 
erators Qi^, Pp, which we define as [12| (note that here 

(Q^) = (Pp) = 0) 

I?« = (RR') - X« = f p)- kM"" , (17) 



M^^ RR^ - (i?i?* 



1 

-1 



(18) 

{Qfj.Qv} : 



where R = {Qx,Qy, Px, Py f and hence Q^i. 
Ptj,. = {Pt. Pu) and L^, = {Q^P, + P.Q^.)/2. 

Eq. (fT7|) is unitaril^elated to the non-negative bosonic 
contraction matrix [l^, [3l| 



V= (ZZ^) -M 

M = ZZ^ - {{Z^fzy = ( -1 



F G _ 
G l + F 

1 



(19) 
(20) 



where Z = {bx,by,bl,blY, F^i, = {blb^), G^i, = {b^b^) 
and 5^ = (Q^ + iP^)/^/2. Since R = UZ, with U = 
we have = UVU^ and = UMU"^ . 
In both the vacuum and the thermal state ([T6| . we 
have 



where f' = in the vacuum state and 

1 d\nZ 1 _ 1 
'^~dX^ " 2 ~ e/S^M - 1 ' 



f 



(21) 



(22) 



in the thermal state. By inverting Eq. ([7]), we then obtain 

(QxQy) = (PxPy) = L^^ = for fi = x,y, and 



{qD 


-<<>+(lA,)^<^->' 


(23) 






(24) 






(25) 



where 



(26) 



case (C) corresponds to k'^ 



h' 



These averages provide all the elements of (IT71) . The 
symplectic eigenvalues of T)^ and T) are coincident and 
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given precisely by and -1 - (Eqs. jUll-dMl)), with 
physical states corresponding to > 0. They are just 

the standard eigenvalues of the matrix DM = Ij^-^p)^ 

or equivalently, V'^M^ = UVMU^ = iCp%) - I f^- 



D. Vacuum entanglement 

The entanglement of the vacuum |0') is a measure of 
its deviation from a product state |0a:) ig) |0y). It can be 
quantified through the entanglement entropy [l3| , which 
is just the von Neumann entropy of the reduced state 
Pfj, = Tr_p p of any of the modes (/x = x,y), since for a 
pure state (p — |0')(0'|) they are isospectral. The state 
Pp, is a gaussian mixed state completely determined by 
the reduced 2x2 covariance matrix 



(27) 



whose symplectic eigenvalues are = det"^/^[X'^ 
^TM^] — ^ and —1 — /p. Here L^^ = and hence, 



(28) 



which is just the deviation of the mode uncertainty from 
its minimum value. The entropy of p^ is, therefore, that 
of a boson system with average occupation /^: 

S{pp) = -Ti-pplogpp = h{U), (29) 
h{f) = -/log / + (1 + /)log(l + /) , (30) 

which is just a positive concave increasing function of 
/p. The vacuum is then entangled iff /^ > 0, with S ~ 
-/^(log - 1) foi' /m ^ and S" « log /^ + 1 for /^ ^ cx) 
(for log = In). 
In the vacuum case {f^ = 0) Eqs. (I23l)-(l28l) lead to 



which is independent of p, where 




(31) 



(32) 



denote the arithmetic and geometric averages of the orig- 



inal oscillator frequencies uj,. 



kfj^. Entanglement is 



thus completely determined by the ratios w/w and Hig/ui 
(with ujg/ijj < 1), or equivalently, uj/uJx and ujyjujx- It 
is then non-zero ^ uj > if Ux ^yi i-^-i w 7^ 
(anisotropic case). In sector A, the are positive, 
whereas in B they are both imaginary, implying 




(fcp < 0) . 



(33) 



E. Thermal entanglement 

For a mixed bipartite state, like the thermal state ([16]) 
at T > 0, entanglement is a measure of its deviation 
froma.epa™6Zestate[4i,i.e.,from a convex combina- 
tion of product states ps = QaPx ^ Py: where qa > 0, 
la — ^- Such states can be created by local op- 
erations and classical communication. For a two-mode 
gaussian mixed state, entanglement can be quantified by 
the negativity fl3\ , which is minus the sum of the nega- 
tive eigenvalues of the partial transpose p*" of the total 
density matrix p, measuring then the degree of violation 
of the PPT criterion [l^ [13 by the entangled state. For 
a two-mode gaussian state, a positive negativity is a nec- 
essary and sufficient condition for entanglement Q . 

Partial transposition with respect to y implies the re- 
placement Py — >■ —Py in the full covariance matrix p7)) 
[sIItI-Q, leading to a matrix V^. The negativity can then 
be evaluated in terms of the negative symplectic eigen- 
values of this matrix, which will have eigenvalues /^ and 
-1 - /^, with /^ > -1/2 [ll|. Replacing L^y by -L^y 
in p7|) . we obtain here 



/± = V"±V«2-/3'-^ 



^i(a + /3)±yi( 



(34) 



where a and /3 can be expressed in terms of the local 
symplectic eigenvalues /^ ([25]) and the global ones /^ 



" = K(Q'>(^'> + (QDiPy)) + {Q.Py){QyP.) 



(35) 



= ^JiiQ.Py)' ~ {Ql){P^mQyPx? - (QDiP^)) 

= n,(/; + ^)- (36) 

(Note that if a is replaced by a = ^{{Ql){Px) + 
{Ql){Py))-{QxPy){QyPx), Eq. m becomes /4). While 
ffj. depends on p for T > 0, Eq. (|35p depends just on the 
sum 



Eu(/m + hf - 2O-^(i0'i+O^) (1 + 



^[0 



.2 ,-.2\ 



The negativity can then be expressed as [l^ 
N = ]-{Ti\p'y\ - 1) = ^[TT ^ - 1] 



Maxf' 



-/- 



1 + 2/- 



^,0], 



(37) 



since only /_ can be negative. The entanglement condi- 
tion /- < leads to a> ^ + 2/3^ or 

E + /m) > e /m(i + /;) + 2 n /m(i + /;) - (38) 
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which imposes a temperature dependent lower bound on 
the average local occupation. 

In the vacuum case = 0, Eq. ([55t implies just > 
0, while Eqs. (I311)-(I371) reduce to 



f-^U- ^JUU + 1) , N = f^ + .^/^(/^ + 1) , (39) 

with given by Eq. pip , in agreement with the general 
results for pure gaussian states 12]. Both — /_ and N 
are again concave increasing functions of at T = 
and can be taken as alternative vacuum entanglement 
measures. 



F. Quantum Discord 

Quantum discord [l^ is essentially a measure of the 
deviation of a bipartite mixed quantum state from a clas- 
sically correlated state, i.e., a state diagonal in a standard 
or conditional product basis. For a general bipartite sys- 
tem A+B, the quantum discord can be defined as the 
minimum difference between the conditional von Neu- 
mann entropy of A after an unread local measurement 
Mb in B and the original quantum conditional entropy 
SiA\B) = S{A,B) - S{B) 13: 

= MAnY,PjS{pA/j) - [S{pab) - S{pb)] , (40) 

Mb "^"^ 



where, for a measurement Mb based on local projectors 
Pj {J2j P] = Ib), Pj = T^PabIa "X) Pj is the probability 
of outcome j and p^/j the reduced state of A after such 
outcome. Eq. PO]) can be also expressed as the mini- 
mum difference between the original mutual information 
I{A : B) = S{A) — S{A\B), which measures all correla- 
tions between A and B, and that after the unread local 
measurement, S{A) — J2j Pj'^iPA/j), which contains the 
"classical" part of the quantum correlations [l^, [2^ . 

For a pure state {p\b — Pab) both S{pab) and 
S{pA/j) vanish and reduces to the entanglement en- 
tropy S{pb) = S{pa), with = For a 
mixed state, however, is not an entanglement mea- 
sure, being in fact non-zero for most separable states [l^l 
and vanishing just for those separable states of the form 
Pc — J2jPjPA/j ^ Pj (classically correlated with respect 
to B), which remain unaltered after the local measure- 
ment AIb- In general, ^ for mixed states. Hence, 
for a bipartite system with a non-degenerate ground state 
in a thermal equilibrium state, like the system under 
study, differences between quantum discord and entan- 
glement, and between and , will arise only at 
finite temperature. 

The exact evaluation of involves a difficult mini- 
mization over all local measurements Mb- Nonetheless, 
for a two-mode gaussian state, a minimization restricted 
to gaussian measurements was recently shown to be an- 
alytically feasible [21I, [g^, 113 • -^'o'" such measurements in 



the present system, Eq. (|40)) becomes, choosing B = y 
and using Eqs. (gl]), ((29)) - (15(11) . 

Dy = MiuM^/f ") - [/i(/;) + hif'_) - h{fy)], (41) 



where /; 



My 



B 



^ 4(/, + 



denotes the symplectic eigenvalue of the co- 
variance matrix I?*^" associated with Px/j, which de- 
pends on the 2x2 covariance matrix determin- 
ing the local gaussian measurement My [2ll, [12] ■ The 
final result was provided in [2^ and can be fully ex- 
pressed in terms of the local invariants A = ^{fx + \Y , 
,C^-2E,(/; + i)^-(/, + i)^and 
D = n.. 4( f,', -I- i)^, which determine the quantity i^Min 

of [221 , with M-VCLMyfx ^^ — ^V-E-Min — \- It can be shown 
that if > 1 the two-mode gaussian state is entangled 
[21I [22} . Moreover, the only two- mode gaussian states 
with D'y = are product states [22|. The expression for 
(local measurement in x) is obviously similar (x ^ y 
in previous formulas). 



III. RESULTS 
A. Vacuum entanglement 

Let us now analyze the main features of Eq. ([HT]) . We 
first consider fixed in ([4]) (charged particle in a mag- 
netic field). In the isotropic case uJx = i^y, ijJ = u)g and 

= V w. There is no entanglement since Lz commutes 
in this case with H and leaves the isotropic product vac- 
uum invariant. For Icj^ - 



16^2(^2+^2) 



< w, Eq. ([311) leads to 



iiJx - LOyf + OiiiOx - LOyf) , (42) 



indicating a quadratic vanishing of in this limit. En- 
tanglement also vanishes for uj ^ (no coupling), where 



, ^ 1 n 1 \ , ,2 



7 (^--)-^ + 0(c.-). (43) 



On the other hand, for uj — >■ 00, a remarkable feature 
is that approaches a finite limit, which depends just 
on the anisotropy ujy/uJx'- For a; ^ |<D| Eq. (|31l) leads to 



fu 



— - 1 



(44) 



In sector A, f^ and hence S{pfj) are then increasing 
functions of lj (Fig. [2|) . Mode entanglement is then en- 
hanced just by increasing the field, although it will sat- 
urate for strong fields. This saturation is a consequence 
of the balance between the oscillator part and the cou- 
pling ujLz in dD), as /c^ = k^^ -\- uP' also becomes large, 
reducing (Q^): For w 00, (Q^) ~ -^'^'^ ~^ while 
(p2) ~ 2^0; 00, leading to the finite limit (|44)) . 

In contrast, and hence entanglement, will diverge at 
the edges of the dynamically stable region. For instance. 
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FIG. 2. Entanglement entropy S = S{p^) (Eq. ((29])) between 
the two modes as a function of u in the vacuum of Hamilto- 
nian ([TJ at fixed fc^ = oj^, for ratios ujy/ujx = 0.1, ... , 0.5 and 
<^o = \i^x\- The top panel corresponds to sector A (attractive 
potential > 0), the bottom panel to sector B (repulsive 
potential < 0, |a;| > Wc), where S diverges for ~> cjc 
(Eq. In both cases 5* approaches the same finite limit 

for uj — >■ oo. 



if — ^ 0, UJ whereas ojg — )■ 0, implying cx 




1.5 



_1 1 1 1 1 1 1 


1 1 II 1 II 1 1 1 


/-l.l 
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FIG. 3. Entanglement entropy between the two modes as 
a function of cj in the vacuum of Hamiltonian ([T]) at fixed 
k'ij_ — ijj'u^ , for ratios uj'y/uj'^ — 0.1, 0.3, . . . , 0.9 (top) and 
{Lo'y/uj'^Y = -1.3, -1.2, . . . , -0.7 (bottom), with uj'q = w^. 
The top panel corresponds to fc^ > 0, with the positive defi- 
nite sector A on the left, where S diverges for uj — >■ oj^j = uj'y, 
and the non-positive sector Bi on the right (a; > = ^x), 
where S diverges at cj^ and vanishes for cj — > oo. The 
bottom panel corresponds to fc^ > and fc^ < (sector 
B2), where S diverges for lj — )■ a;^ and also lj — )■ oj^s (Eq. 



(O) if {u}y/uj'^) < —1. In the critical case (iu'y/uj'^y = — 1 
{k'y = —k'^), and S saturate for large uj (Eq. (fS^ ). 



- 1 



(45) 



and hence S{p^) ~ ^\og{uJx/iOy) plus constant terms. 
This divergence stems from that of (Qy) (or (P^)) in 
this limit, with (Py) and (Q^) remaining constant (Eqs. 

In the repulsive sector B, diverges for a; — > \uj\ — Uc 
(Eq. in])), where both (Q^) and (P^) diverge: 




(46) 



It is then seen that here and hence S{p^) decrease 
as Lo increases from uj^ (Fig. [21 bottom panel), i.e., as 
the system becomes dynamicaUy stabilized by the field, 
reaching for w — 00 the same previous limit (l44l) . At 
fixed \ky/kx\, the vacuum entanglement is then strictly 
larger in the unstable sector B (fc^ < 0). 

At fixed k'^ (rotating potential) the behavior with fre- 
quency is quite different (Fig.[3l). We should now replace 



= Jk 



in Eqs. (l3T|)-(l32l). For 
entanglement. For ~ ujy \ ^ ujx, we have 



cjy there is of course no 



(c.;-c.;)2 + o((c.:^-c.;)4). (48) 



0, where Eq. (|43 



Entanglement also vanishes for uj 
still holds {u'^ — utf^ at w = 0). 

On the other hand, as uj increases, increases rapidly 
and in contrast with the previous case, it diverges for 
UJ — > uj'^i (Eq. (fT2)) ). where, assuming uj'^^ = ( ■' 



ujy < UJ^ 



fu 



uj^iuj'i-uj't^) 



2(c.;-^)(3<^+0 



2^2 



1 



(49) 



implying 5(/9^) w j \n[uj'y / {uj'y — uj)] plus constant terms. 
In this limit (Q^) and (P^) diverge while (Ql) and (P^) 

As UJ increases 



0. 



stay constant, as ujy = 
further, the system enters the instability window, al- 
though for UJ > uj'c2 (Eq. ([E])), it recovers a discrete 
spectrum, entering sector Bi. For uj — >■ w^ji /m diverges 



(47) as in dUl), with uj' O uj'^ if 



c2 



^ UJ^. 
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In sector Bi, and hence the entanglement decrease 
as uj mcreases, vanishing for cj —> cxi, in contrast with the 
behavior at fixed in sector A. In this Umit the vacuum 
of H becomes now that associated with uL^, which is an 
isotropic product gaussian state with Lz ^ 0, and hence 
zero entanglement. (Q^) and (P^) stay then finite and 
their product approaches minimum uncertainty, leading 
to 



(50) 



In the unstable domain ,62 , the behavior with uj is the 
same as in Bi when fc^ > and — fc^ < fc^ < 0. However, 
for fc^ > and — 3fc^ < k'y < —k'^, we also have the upper 
instability limit P3)) {uj'^^). In this case first decreases 
with increasing w, reaching a minimum, but then starts 
again to increase, diverging for uj — > where now both 
{QD and (P^) diverge, leading to 



^c3 



- I 



(51) 



We then obtain different 0(aj— Wc)^^^"' and 0(a;c— w)^^/^ 
divergences of at the stability borders uj'^ and cj^g re- 
spectively. 

In the special critical case k'y = —k'^ [uj'y — iuj'^), where 



^c3 



00, Eq. (jSTj) leads to 



1], 



(52) 



and hence to a finite asymptotic limit — ^(v^ — 1) 
for a; — > 00, in contrast with (I50p . as also appreciated in 



Fig. El In this limit (Q^) diverges whereas (P^) vanishes, 



the product approaching 1/2. Hence, as uj increases, 
vanishes if < —k'^, saturates if fc(^ = —k'^, and diverges 
(at UJ = w^g) if -3fc; < A:; < -fc^. 

The behavior of (and hence S{p^)) with w is quali- 
tatively similar to that of the average angular momen- 
tum (iz)- At fixed kf^, the latter also saturates for 
cj — > 00 {{Lz) — >■ {uj/ujgY — 1) and diverges for w — > 
((L^) cx {uj — Wc)~^^^), whereas at fixed k'^ it diverges for 
UJ cj^j ((-^2) oc (w — w^)"-^/^ for z = 1, 2) and vanishes 
for cj — >■ 00. Entanglement is then an increasing function 
of \{Lz)\ at fixed k^ 



or A:;,^ 



as seen in Fig. SI although 
it is not fully determined by K^z)], as the latter is not 
invariant under local transformations (in contrast with 
/p). At fixed {Lz), higher ratios ky/k^ < 1 originate 
a higher entanglement (Fig. For small uj, (Lz) oc uj 
and hence, oc (i^)^ for small (Lz) in sector A. How- 
ever, at fixed k'^, (Lz) also vanishes for large w, where 
{Lz) oc uj~^. Hence, in sector Bi and according to Eq. 
(ISUl) . oc {Lz)^^^ for small (iz), leading to an infinite 
initial slope (dotted line in Fig.^]). At fixed (Lz) and k'^, 
entanglement is then stronger in the unstable sector Bi 




FIG. 4. The entanglement entropy between the two modes as 
a function of the scaled average angular momentum (Lz) = 
{Cz)/h (Eq. (0)), in the vacuum of ^ at fixed fc^ — ujf^ (solid 



and dashed lines) and at fixed fc^ 



(dashed-dotted and 



dotted lines), for selected ratios ujy/ujx and uiy/u'^ and vari- 
able UJ. The solid and dashed-dotted lines (in red) correspond 
to sector A (fc^ > or fc^ > and \uj\ < uj'y), the dashed and 
dotted lines (in blue) to sector B (fc^ < 0, \uj\ > cjc) and Bi 
{k'l^ > 0, \uj\ > uj'x). S is in all sectors an increasing function 
of |{L.)|. 



{uj > uj'^2)- An exceptional behavior occurs in the critical 
case k'y — —k'^ (Eq. ([5^ ). where for w — >• 00, (Lz) oc 
vanishes while remains finite. In this special limit 
there is finite entanglement with vanishing angular mo- 
mentum. On the other hand, close to the divergences, 
(X (Lz) {uj ^ UJ,) or (Lz)i/2 (a; ^ i = 1,2), 
implying S{p^) oc In (Lz) for large (Lz). 



B. Thermal entanglement 

Let us now examine the thermal entanglement in the 
stable sector A. We first depict in Fig. [5] the limit tem- 
perature for entanglement Te, determined from the con- 
dition /_ = (equality in Eq. I\'68\i ). This temperature 
remains finite for all values of k^ or k'^, including the 
edge of the sector (A:^ — >• or |w| — > y^), where the 
vacuum diverges. At the edge, A_ — > and hence 
a finite T already gives rise to a spread over all energy 
levels (/L — 00), which diminishes and eventually kills 
the entanglement. A related fundamental effect is that at 
finite T > 0, entanglement does not diverge at the edge, 
but stays finite or vanishes, depending on the value of T. 

More precisely, for ujy and fixed uj^ > 0, 
A+ — > Auj"^ + ujI whereas A_ 



« r/A_ « TX+/{uj^ujy). 
lead to 



ujyUJx/X+, implying 
Hence, in this limit Eqs. 



TA* (1 



V'+)' 



2^\/c.2(2^2A^(H-2/;)+rc.2) 



- 11 



(53) 



which remains finite and above —1/2 if T > 0. This 
implies a finite negativity in this limit if T > 0, with 
N oc for T ^- according to Eq. Therefore, 

at finite temperature the vacuum divergences of the en- 
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FIG. 5. Top: Scaled limit temperatures for entanglement Te 
at fixed fc^ (top) and at fixed k'^^ (bottom), as a function of 
frequency for different ratios ky/kx = 0.01, 0.1, 0.9, 0.99 
(top) and k'y/k'x — 0.05, 0.1, 0.95 (bottom), with ujo = 
ijJx = oj'o = iu'x — \fkJ^. At fixed fc,t, Te vanishes for 

large cj (Eq. (|55|l ') whereas at fixed fc^, it approaches a finite 
value at the upper stability limit = \f^)- 




FIG. 6. Top: Negativity as a function of frequency at in- 
creasing temperatures T/ljo = 0.001, 0.03, 0.06, 0.15, at 
fixed kp, (top), for ky/kx ~ 0.2, and at fixed fc^ (bottom), for 
k'y/k'x = 0.5. At T > 0, A'^ is finite, and non-zero just within 
a finite frequency window. 



tanglement can be only probed indirectly, through the 
behavior of near the edge at sufficiently low T. 
In addition, Eq. (j53p entails a finite limit temperature 
Te, obtained from the condition /_ = in ([5^ : 



Te 



2(1 + 2/;)^2^gA+ 
(l + 2/;)2A4 



(54) 



which is a transcendental equation for Te (/+ depends 
on Te). The maximum limit temperature T^ at fixed 
kx = oj^ or T'e at fixed k'^ — oj'J^ , is in fact obtained in 
this hmit (ujy = or ujy = lo): At fixed k^, ~ 0-24:UJx, 
attained at w « 0.38uJx, while at fixed k'^, T''^ « 0.23cj^, 
attained at cj « 0.28uj'x- 

At fixed fc^, the limit temperature Te as a function of 
uj exhibits first a maximum and then vanishes for w — > cx) 
(top panel in Fig.[S]), i.e., in the limit where the vacuum 
entanglement saturates. The reason is that A_ also van- 
ishes for a; — c» (A_ ~ ujxt^y / i'^^)) , implying Te c>c uj^^ 
in this limit: 



Te 



2uj In 



(55) 



In fact, for a; — )■ c» and fixed fc^ (with ky < kx), Eqs. 
dM!)-® lead to 




such that /_ = for 1 + 2f'_ = ujx/<^y, which leads to 
Eq. ([55t . On the other hand, for a; — > 0, T vanishes 
logarithmically {Te oc — 1/lnw) and the same occurs for 
ujy — >■ uJx, since in these limits A± remain both finite 
whereas the negativity vanish. At fixed T, we then obtain 
a finite frequency window for entanglement, which nar- 
rows for increasing temperature or decreasing anisotropy, 
as seen in Fig. 5 and also Fig. [6l where the negativity 
([57)1 is depicted. Let us remark that entanglement ceases 
to be correlated with (L^) as the temperature increases 
{{L,) (X LoT for high T). 

At fixed k'^, the behavior of Te and N look quite 
different, as now lo is bounded above by w'y (assuming 
uo'y < uj'^). For oj oj'y, Te is then determined by Eq. 



((54)) with UJx ^ y oj'x —uS^, and remains finite. Actu- 
ally, as verified in the bottom panel of Fig. [51 Te acquires 
in this border its maximum value as lo increases at fixed 
k'^ if uj'y < uj'y^ « 0.28a;^, while if uj'y > ujy^ the maximum 
is attained at an intermediate frequency. Consequently, 
at fixed T < T'e there is again entanglement within a 
certain frequency window, which extends up to the sta- 
bility border uj = uj'y ii u!y < uj'y^ or T < Te ai uj = uj'y 

(bottom panel in Fig. 16]). The absolute maximum T'^e is 
obtained at this border precisely at uj'y — uj'y^. For w — > 



(56) 

or uj' uj'^, Te decreases again logarithmically. 



9 



0.4 


1 1, 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 


0.3 


■ 1 ~~~~ 




0.2 






0.1 


■ 







■ . . N 


■ 






0.2 0.4 0.6 0.8 


1.0 1.2 



FIG. 7. Negativity (A''), single mode entropy — S{px) and 
the quantum discords and as a function of temperature 



for uj = Loo 



kx and ky 



0.2kx. 



that C cx T-^: 



(57) 



with a similar expression for after replacing ujx by ujy . 
Hence, for high T becomes independent of a;^, with 
> Dy asymptotically if LOy < tOx, as verified in Fig. [T] 
We also remark that the discord remains finite for T > 
in the whole sector A. 



IV. CONCLUSIONS 



C. Comparison with the quantum discord 

We finally compare in Fig. [7] the thermal behavior of 
the negativity with that of the gaussian quantum discord 
(Eq. (HI])) and also D^. For reference we have also 
plotted the entropy of one of the modes (x), no longer 
a measure of entanglement, just to indicate its coinci- 
dence with both Dy and Z?^ for T 0. While at T = 
the negativity is just an increasing function of the en- 
tanglement entropy (Eqs. (p9)) - (p9)) ) and hence of the 
quantum discord, the behavior for T > is quite dif- 
ferent. Although exhibiting a similar initial decreasing 
trend (essentially due to the initial increase of the total 
entropy S{pab) in (001)) tbe gaussian discord starts then 
to increase (due to the increase in the first term of (|40p ). 
vanishing only asymptotically for T — >■ oo. Such revival 
of the discord with increasing T was also observed in 
spin systems 0, [EOl , and reflects the presence of quan- 
tum correlations in the excited eigenstates, which lead 
at these temperatures to a separable yet not classically 

thermal state. Since 
[22i. one can ensure 



correlated (in the sense of sec. IIIF 
D^^ > 1 implies entanglement [21 
here that < 1 after the vanishing of the negativ- 
ity (T > Te), although this may not prevent from 
reaching a higher value than at T = at some interme- 
diate temperature, as seen in Fig. [T] For T — oo we 
actually obtain, from Eq. (01]) and the expression of [13] , 



We have analyzed the entanglement induced by an 
angular momentum coupling on two harmonic modes. 
Full analytic expressions for the vacuum entanglement 
entropy and the thermal negativity were derived. The 
model exhibits a rich phase structure and admits dis- 
tinct physical realizations (particle in a magnetic field 
in an anisotropic harmonic trap, or particle in a rotat- 
ing harmonic trap) , which lead to different entanglement 
behaviors with the relevant control parameter. For in- 
stance, in sector A (stable vacuum), entanglement sat- 
urates for strong fields in the first case, but diverges at 
a finite frequency in the second case. Vacuum entangle- 
ment diverges at the onset of instabilities, being corre- 
lated with the average angular momentum and reaching 
higher values in unstable domains dynamically stabilized 
by the field or rotation. In contrast, thermal entangle- 
ment is finite, and non-zero just below a finite limit tem- 
perature within a reduced frequency window, diverging 
only for T — > at the instability borders. We have also 
shown that after a short initial common trend, the ther- 
mal behavior of the gaussian quantum discord becomes 
substantially different from that of entanglement, van- 
ishing only asymptotically. A deeper investigation of the 
discord and other related measures of quantum correla- 
tions [48, 51] in similar systems is being undertaken. 

The authors acknowledge support from CONICET 
(LR) and CIC (RR) of Argentina. 
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